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Abstract 



In this paper, Neumann cracks in elastic bodies are considered. We establish a 
rigorous asymptotic expansion for the boundary perturbations of the displacement 
(and traction) vectors that are due to the presence of a small elastic linear crack. 
■ The formula reveals that the leading order term is e 2 where e is the length of the 

, crack, and the e 3 -term vanishes. We obtain an asymptotic expansion of the elastic 

<-j ' potential energy as an immediate consequence of the boundary perturbation formula. 

The derivation is based on layer potential techniques. It is expected that the formula 
would lead to very effective direct approaches for locating a collection of small elastic 
cracks and estimating their sizes and orientations. 



Mathematics subject classification (MSC2000): 35B30, 74B05 
Keywords: elastic crack, expansion formula, boundary perturbations 

£q ■ 1 Introduction 

The displacement (or traction) vector can be perturbed due to the presence of a small 
£NJ ■ crack in an elastic medium. The aim of this paper is to derive an asymptotic formula for 

the boundary perturbations of the displacement as the length of the crack tends to zero. 
The focus is on cracks with homogeneous Neumann boundary conditions, i.e., perfectly 
insulating cracks. We consider the linear isotropic elasticity system in two dimensions and 
assume that the crack is a line segment of small size. The derivation of the asymptotic 
formula is based on layer potential techniques. 

The paper extends recent asymptotic results that have been used for an efficient imag- 
ing of small defects. In [8] an electrostatic model, where the crack is perfectly conducting, 
was considered and an asymptotic expansion of the boundary perturbations that are due 
to the presence of a small linear crack was derived. The asymptotic formula leads us 
to efficient algorithms to detect cracks using boundary measurements [3l [5]. Their res- 
olution and stability of the algorithms with respect to medium and measurement noises 
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were investigated in (2j H] . There were some work on boundary perturbation due to the 
presence of small inclusions in linear elasticity; the effect of small inclusions on boundary 
measurements has been studied in jSJE]. The effect of thin elastic inclusions on boundary 
measurements was quantified in |101 1 12 j . Direct reconstruction algorithms for locating 
small or thin elastic defects were developed in [IJ [Til HS1 US]. We emphasize that the 
results of this paper (on cracks) can not be obtained as a limiting case of thin inclusions. 

The results of this paper reveals that the leading order term of the boundary pertur- 
bation is e 2 where e is the length of the crack, and its intensity is given by the traction 
force of the background solution on the crack (see Theorem 14.11 ) We also prove that the 
e 3 -order term vanishes. By integrating the boundary perturbation formula against the 
given traction, we are able to derive an asymptotic expansion for the perturbation of the 
elastic potential energy, which is an improvement over already existing results [26, 22J (see 
the discussion at the end of Section 5). 

The boundary perturbation formula derived in this paper carries information about 
the location, size, and orientation of the crack, and we expect, as in the electrostatic case, 
that the formula will provide a powerful tool to solve the inverse problem of identifying the 
cracks in terms of boundary measurements. The implementation of imaging algorithms 
based on the present expansion and the analysis of their resolution and stability will be 
the subject of a forthcoming paper. 

The paper is organized as follows. In section [21 a representation formula for the 
solution of the problem in the presence of a Neumann crack is derived. Section [3] is 
devoted to making explicit the hyper-singular character involved in the representation 
formula. Using analytical results for the finite Hilbert transform, we derive in section H] 
an asymptotic expansion of the effect of a small Neumann crack on the boundary values 
of the solution. Section [5] aims to derive the topological derivative of the elastic potential 
energy functional. Appendix contains technical calculation of the double layer potential. 

2 A representation formula 

Let VL C M? be an open bounded domain, whose boundary d£l is of class C 1,a for some 
a > 0. We assume that U is a homogeneous isotropic elastic body so that its elasticity 
tensor C = (Cijki) is given by 

Cijki = ^Sij5ki + fi(SikSji + SuSjk) (2.1) 

with the Lame coefficients A and fi satisfying fj, > and A + /x > 0. Let 7 E C O be a small 
straight crack with size e, located at some fixed distance do from d£l, i.e., 

dist(7 e , dQ) > do. 

We denote by a unit normal to j £ . 
Let 

^ := U : Wj + Oji/Ji = 0, 1 < i, j < 2 j, 

or equivalently, 



^ = span 
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Introduce the space 



Ly(d£l) := if 6 L 2 (dU) : f f-ipda = ferall^G^}. 

I Jan J 

Let u e be the displacement vector caused by the traction g £ L\{dVL) applied on the 
boundary d£l in the presence of 7 e . Then u e is the solution to 



V-a(u e ) = mfl\%, 

a(u e ) n = g on <9f2, (2.2) 
a(u E )e ± = Q on7 e , 

where n is the outward unit normal to <9f2 and c(u e ) is the stress defined by 

a(u £ ) = CV s u E := ^C(Vu E + VuJ). (2.3) 

Here, V s u £ = ^(Vu e + VuJ) is the strain tensor and the superscript T denotes the 
transpose of a matrix. Note that the functions in \& are solutions to the homogeneous 
problem (|2.2h with g = 0. So we impose orthogonality condition on u e to guarantee the 
uniqueness of a solution to (|2.2|) : 



/ u £ • ip da = for all (2.4) 
Jdn 

Let uo be the solution in the absence of the crack, i.e., the solution to 

fv • <7(uo) = in O, 
lcr(uo)n = g on <9f2, 

with the orthogonality condition: uoIqq G L|,(<9f2) (or equivalently, (|2.4p with u e replaced 
with uo). 

It is well-known that the solution u £ to (|2.2p belongs to H 1 ^ \ j E ). In fact, we have 

M^nyy.) < C ( 2 -6) 
for some C independent of e. To see this, we introduce the potential energy functional 

J £ [u\:=-\f a(u):V s u. (2.7) 

The solution u e of (|2.2p is the maximizer of J e , i.e., 

J £ [u £ ] = max J e [u], (2.8) 

where the maximum is taken over all u G H 1 ^ \ -y e ) satisfying c(u) n = g on dQ and 
cj(u) e -1 - = on 7 e . Let v be a smooth function with a compact support in Q such that 
cr(v) e 1 - = — ct(uo) e 1 - on 7 e . We may choose v so that J e [v] is independent of e. Since 
> J e [u e ] > J e [uo + v], we have 

l|V s u £ || i 2( Q \ 7e ) < C. 
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We then have from the Korn's inequality that there is a constant C independent of e such 
that 

||U £ - Uo||h1(0\7«) < C(l|V S (u £ - Uo)||ra(n\ 7 .) + \\ u e - u o|li/i/2(9C))- ( 2 - 9 ) 

Since ||u £ — ^■o\\H 1 / 2 (dn) is bounded regardless of e as we shall show later (Theorem 14. ip . 
we obtain (12.61) . 
Let 

<^ e (x) := u £ | + (x) - u £ |_(x), xE7 £ , (2.10) 

where + (resp. — ) indicates the limit on the crack 7 e from the given normal direction 
(resp. opposite direction), i.e., 

u±(x) : = lim u(x ± te^). 

i->0 

We sometimes denote c(u) n, the traction on dil (or on j £ ), by dn/dv, i.e., 

:= A(V • u)n+ u(Vu + Vu T )n on dtt. (2.11) 
If $ = (3 , i :; )2x2 is the Kelvin matrix of the fundamental solutions of Lame system, i.e., 
^(x):=^<%log|x|- Jj-^, x/OGK 2 , (2.12) 



where 

. A + 3/j, . „ A + /i . . 

yl = — ; — — r- and £ = — - — — — , (2.13) 

2/x(A + 2^) 2/i(A + 2)Lt)' V ; 



then the solution u £ to (|2.2p is represented as 

u e (x) = / |i(x- y)u £ (y) da(y) - / <f(x- y)^(y) dtr(y) 
Jan o'^y Jan 

-/ ^— (x-y) Ve (y)dff(y), xe^\ 7e . (2.14) 



'7s ^"y 

The solution uq to (|2.5p is represented as 



, , f d& . . . . . . f , , 5uo . . . . 

u °( x ) = / ^— ( x -y) u o(y)^(y) - / $( x -y)i^-(y)*Ky)- 

Jan o^y Jan 

Let 

w e := u e - u . (2.15) 
Since = on <9f2, by subtracting above two identities, we have 

w £ (x) - X> n [w fi ](x) = -De^Kx), xe!l, (2.16) 

where the double layer potentials Vq and P e are defined by 

2?n[w e ](x) := / — (x-y)w e (y)cHy), x ^ O 2 - 1 ?) 
Jan "^y 
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and 



T> e [<p e ](x) := f |^(x-yV £ (y)da(y), x£!l. (2.18) 
Let n x denote the outward normal to <90 at x G d£l and let 

%ML( X ) - lim %M(x + tn x ). 

1 t-s>0+ 

Then, it is well known (see, for example, |6|) that 

%M|_(x)= U J + ^n) M(x), a.e. x G <9f2, 
where /Cq is the boundary integral operator defined by 

/C^M(x) := p.v. / — - (x-y)w £ (y)dcr(y) ! x G dQ, 

and / is the identity operator. Here, p. v. stands for the Cauchy principal value. It then 
follows from ()2. 16j) that 

1 



2 -J + /C n J [w e ](x) = V e [<p e ](x), x G an. (2.19) 
Since — \l + /Cq is invertible on L|,(<9ft) (see, for instance, [13]), we have 

w £ (x)=/ A(_lj + x: n )- 1 [$(.-y)]( X )^(y)d,7(y), XE90. (2.20) 
Note that 

(-i/ + /C Q )- 1 [<D(.- y )](x) = N(x, y ), xGffl, yGft, (2.21) 

modulo a function in ^, where N(x, y) is the Neumann function for the Lame system on 
O, namely, for y G ft, N(x, y) is the solution to 

r V-a(N(-,y)) = -S y I in ft, 

(2 22) 

ff (N(-,y))n = -— I on 3ft, 

subject to the orthogonality condition: 

/ N(x,y) • *0(x) da(x) = for all if ef. 

Here, I is the 2 identity matrix. See [6l [T7] for properties of the Neumann function and 
a proof of (pT2"T|) . Thus we obtain from (^0|) that 

u £ (x) = u (x) + / — N(x,y)</> £ (y)da(y), xeffi. (2.23) 



We now describe the scheme to derive an asymptotic expansion of u e — uo on dft. Since 



Qrf- = cr(u e )e- L = on 7 £! we use (12. 16|) to obtain 



^ + ^-V n [^ £ ] = ^-V £ [ V£ ] on 7e . (2.24) 
ov du ov 

We solve this integral equation for (p £ and then substitute it into (|2.23p to derive an 
asymptotic expansion of u £ as e — > 0. 
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3 Derivation of an explicit integral equation 

In view of (I2.24|) . we need to compute -J^- (j^-fa — y)^ on j 6 . As before e 1 - is the unit 

normal to j £ , and denoted by e 1 - = (ni,n 2 ). It is worth mentioning that n\ and n 2 are 
constant since "y £ is a line segment. It is convenient to use the following expression of the 
conormal derivative: 

(3.1) 



7T = r 9 u > 
av 



where the operator T(d) = T(di, d 2 ), where dj = is defined by 



(A + 2/x)niCi + /xn 2 6 + Arii£ 2 

An 2 6 + (Jtnib /in^i + (A + 2^)n 2 £ 2 



(3.2) 



We first obtain the following formula whose derivation will be given in Appendix lAl 
For x/y, we have 



( 



\dv. 



1.) 



aSa + b 



{xi - yi)(xj - yj) 



where 



E 

1=1 

nj(xj - yj) - m(xj - yj) 

I* - y| 



ni(xi - yi) 



27r(A + 2/i)' vr(A + 2/x) 

Let Vij = (J^-(x — y))ij for convenience and let 



d { d§ 



W{*-y) : — — ix 



V 



(3.3) 
(3.4) 

(3.5) 



Then one can use (13.21) to derive 



W(x-y) n 
W(x - y)i 2 
W(x - y) 21 
W(x - y) 22 



ni[(A + 2fi)div u + A5 2 u 2 i] + n 2 {nd 2 vu + fidiv 2 i), 
m[(A + 2fi)d!V 12 + Ad^u^] + n 2 (fid 2 vn + ^idiv 22 ). 
ni(fj.d 2 vn + iidiv 2l ) + n 2 [A9it;ii + (A + 2/i)<9 2 v 2 i], 
ni(nd 2 v\ 2 + fJ,div 22 ) + n 2 [A9iVi 2 + (A + 2/i)<9 2 u 22 ]. 



Since the crack which we consider is a line segment with length e in the domain Q C 
we may assume, after rotation and translation if necessary, that it is given by 



7£ = {(xi,0) :-£/2< Xl <e/2}. 



(3.6) 
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In this case, one can check that 

Oi -yi) 2 



d 2 v u 
d\V2i = a 
div 12 = - 

2 V 22 



a + b- 



a + b 



x - y| 2 J |x - y| 2 (x x - yif 

1 o x i - vi x i - vi 



x-y| 
1 



|x - y| 2 \xi - yi| 2 
xi - yi xi - yi 



l x -y| 



(xi - yi) 2 ' 
a 

|x - y| 2 \xi - yi\ 2 ) fa-y^ 2 ' 



(xi - yi) 2 
d\vu = d 2 V2i = d 2 vi2 = div 22 = 0. 



Since e 1 - = (ni,n 2 ) = (0, 1), we have 



W(x-y) 



ii 



n{a + b)- 



1 



1 



\i a 



lib 



W(x-y) 12 
W(x-y) 21 
W(x- y) 22 



(xi - yi) 2 
Hd 2 vi2 + ^div 2 2 = 0, 
Adiun + (A + 2/i)d 2 v 2 i = 0, 
\d\vi2 + (A + 2/j)<9 2 t> 22 

1 1 



Oi - 2/i) 2 (xi - yif 



Xa 



(xi - y\) 2 



+ (A + 2fi)a 



__ 2(A + jj)a 
(xi-yi) 2 (a^i-yi) 2 ' 



that is, 



Note that 



W(x - y) 



7r(A+2^t) 





(xi - yi) 2 

m(a + m) s 





tt(A+2^) 



A + 2^ 4 

where is the Young's modulus in two dimensions. So, we have 

E 1 



W{x - y) 



4-7T (xi - yi) 2 

So far we have shown that if 7 e is given by (13. 61) . then 



I. 



dy, x = (x,0), -e/2<x<e/2. 



(3.7) 
(3.8) 

(3.9) 
(3.10) 



Here the integral is hyper-singular and should be understood as a finite part in the sense 
of Hadamard, which will be defined in the next section. So the integral equation (12.24}) 
becomes 



where 



f(x) = ^°(x,0) + | ; %[w £ ](x,0). 



(3.12) 
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Define 

f £ (x) := f (|x), (3.13) 

and 

A(x) := -<P*{\x), -K x < 1. (3.14) 
Then the scaled integral equation is 

i I My) dy = _i < x < i (3.i5) 

which we solve in the next section. 

4 Asymptotic expansion 

The integral in (|3.15|) is understood as a finite-part in the sense of Hadamard [19, 20]: for 
i/> GC 1 ' Q (-1,1) (0 < a < 1) 



1 Hv) J v r f x ~ s Wv) j , f 1 J 2^(x) 



■ = lim 



l - V) 2 S->o L y_ x (x - y) 2 7^+5 (x - y) 



/ t v? d y 



(4.1) 



Define 

>t^](x):=-/ ,^ y \ 2 dy, \x\<l. (4.2) 

W-i - yr 

It is known (OH]) that 

AMx) = ~nmx), (4.3) 
where % is the (finite) Hilbert transform, i.e., 

nmx)=p.v.- f 1 ^dy. (4.4) 
vr J -i x - y 

More properties of finite-part integrals and principal-value integrals can be found in \17\ 

na [2oi eh [la [27] . 

If V'(-l) = V'(l) = 0, we have from P~3j) that 

.4[<i/>](x) = -W[^(x). (4.5) 

Thus we can invert the operator A using the properties of T~L. The set y, given by 

i 



y=Vp:l Vl-x 2 \p(x)\ 2 dx <+oc\, (4.6) 
is a Hilbert space with the norm 

IMIy= (J^VT^Wxtfdxy 2 



s 



It is well known (see, for example, section 5.2]) that T~L maps y onto itself and its null 
space is the one dimensional space generated by l/Vl — x 2 . Therefore, if we define 



X 



i> € C° ( [-1,1] ) : € y, V(-l) = V-(l) = , 



(4.7) 



where tp' is the distributional derivative of tp, then A : X — > y is invertible. We note that 
X is a Banach space with the norm 



\x = M\L~ + W\\y- 
Using the Hilbert inversion formula (see, for example, section 5.2]), we can check 



(4.8) 
(4.9) 



The equation (|3.15|) can be written as 

4 



A[ip E ](x 
The Taylor expansion yields 



E 



f e (s), -1< x < 1. 



(4.10) 



du s n x du . exd 2 u 



dv y 2 ' dv w 2 <9t^ 
where <9/<9i denotes the tangential derivative on %. The remainder term e\ satisfies 

\ei{x)\ < Ce 2 \x\ 2 , 

and in particular, 



(4.11) 



|eib < Ce\ 



(4.12) 



On the other hand, since 



d 2 N „e 



d e . 

— V n w £ \{-x, 0) v 
di/ 2 ./an ov x dv y 2 



((-x,0),y)w e (y)^(y), 



and 7 £ is away from <9f2, one can see that 



< C||w e || io 



(an)- 



Therefore, we have 



where e satisfies 



M»)-^<o) + f £>) + .(«), 

az/ 2 ata^ 

||e||y < C(e 2 + ||w e || i0 o (9f2) ) . 



(4.13) 

(4.14) 
(4.15) 
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We now obtain from (|4.10|) that 



^(O)^ 1 [!](*) + i^-^A-'Mx) +A- 1 [e)(x) 



(4.16) 



Note that £\{x) = A 1 [e\{x) satisfies 



||£llUr<C|N|y<C U 2 + ||w £ || L 



'(da) 

and in particular, 

ll^i (-1,1) < C ^e 2 + ||w £ || LO o( ar2 )^ . 
It then follows from (USD and (El that 

ip £ (x) 



E 



^(0)71^ + ||g(0)-^5 + £l{x) 



(4.17) 



(4.18) 



Thus we have from (|3.14p that 



du 



ld 2 u 



(0)Ve 2 - 4x 2 + -^-f(0)xVe 2 - Ax 2 + S(i 
ov 2 atav 



(x,0)e 7e , (4.19) 



where £{x) = |£i(~x) satisfies 



||£||z,°o( 7e ) < Cs\ e 2 + ||w e ||x,oo(an) 



(4.20) 



Substituting ([Ug]) into (pT23l) we obtain 

w B (x) = |y A N(X)( y )0 ))7^T^^^ ( o) 

+ 1^ ^N(x,(y,0))£(y)dy.= I + II + III, x G 90. 



(4.21) 



Since 



<9f. 



( N(x, (y, 0)) = ^ 7 N(x, 0) + — N(x, 0)y + 0(y 2 ), 



we have 



/ ^-N(x, (y, 0)) \A 2 - V ^ = t^-N(x, 0) / y/e 2 -4y 2 dy 



+ 



dt 

ire 2 d 



2 9ia 



a 1 r 

^— N(x,0) / yy/e 2 -Ay 2 dy + 0(e 4 
N(x,0) + O(e 4 ), 
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and hence 

'=^4 n(x ' 0) ^ (0)+o(e4) ' ^ 

Here and throughout this paper, 0(e 4 ) is in the sense of the uniform norm on dQ. Similarly, 
one can show that 

U = 0(e 4 ). (4.23) 

So we obtain that 

TT£ 2 d 1VT/ n ^duo fnS , 4 . 



w e (x) = — — N(x, O)-^(O) +0(e 4 ) + ///. (4.24) 
In particular, we have 

||w 6 |Uoo (an) < C(e 2 + \III\). 
But, because of (|4.20p . we arrive at 

\III\ < Ce\\S \\ L oo {le) < Ce 2 (e 2 + ||w £ || ioo(an) ), 

and hence 

II w £||l°°(9C) < Ce (1 + \\we\\L™(dn))- 
So, if e is small enough, then 

l|w £ || L oc (9n ) < C£ 2 . (4.25) 

It then follows from (I4.24j) that 

w £ (x) = ^A N(x , o)^(0) + 0(e 4 ). (4.26) 

We obtain the following theorem. 

Theorem 4.1. Suppose that 7 £ is a linear crack of size £ and z is the center ofj E . Then 
the solution to h2.2\j has the following asymptotic expansion: 

, . . TT£ 2 <9N , , Chin , , <, . 

(u e - u )(x) = ( x ,y) — ° (z) + 0(e 4 ) (4.27) 

Hi OVy y=z OV 

uniformly on x € <9fL Here E is the Young's modulus. 

It is worth emphasizing that in (|4,27p the error is 0(e ) and the e 3 -term vanishes. 
One can see from the derivation of (|4.27p that the e 3 -term vanishes because 7 e is a line 
segment. If it is a curve, then we expect that the e 3 -term does not vanish. We also 
emphasize that (|4.27p is a point-wise asymptotic formula, and it can be used to design 
algorithms to reconstruct cracks from boundary measurements. We can also integrate this 
formula against the traction g to obtain the asymptotic formula for the perturbation of 
the elastic energy as we do in the next section. 

Similarly, if we consider the Dirichlet problem 



V • tr(u £ ) = in \ %, 
u £ = f on dn, (4.28) 

k cr(u £ )e- L = on7 £ 



and denote the Green function of Lame system in by G, then we get the following 
asymptotic expansion of its solution u e . 
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Theorem 4.2. Suppose that j £ is a linear crack of size e, located at z. Then the solution 
to ^4-28ty has the following asymptotic expansion: 



d . „ , tie 2 d 2 G 

— (u e - u )(x) = — — 

ov E ov„ov. 



(x,y) 



x^"y 



5uo 

y=z 9f 



(z) + 0(e 4 ) 



(4.29) 



uniformly on x € <9S7. 



5 Topological derivative of the potential energy 

The elastic potential energy functional of the cracked body is given by (|2.7p . while without 
the crack the energy functional is given by 



J[u ] = -- / cr(u ) : V s u . 
1 Jn 

By the divergence theorem we have 

Je[ u s] ~ J[ u o] = —x [ (u e - u ) • gala. 
1 JdQ 

Thus we obtain from (|4.27p 

J e [u e ]-J[u ] = -^^(z)A / N(x,y)g(x)da(x) + 0(e 4 



(5.1) 



(5.2) 



Since 



we have 



2£ <9zv' 

u (y)= / N(x,y)g(x)d<r(x), y G fi, 

- 2 1 <9u 2 



■4[uJ - ^[uo] 



7T£ 

2£ 



+ 0( £ 4 ). 



(5.3) 



We may write (|5.3p in terms of the stress intensity factors. The (normalized) stress 
intensity factors Kj and ifjj are defined by 



So, we have 
and hence 



Ki(uo,e) := cr(uo)e • e 1 " and AT//(uo,e) := er(uo)e • e. 
cr(u )e- L = Kie 1 - + K n e, 



i>U0 (z) 2 = \a(u )e^ = K 2 + K] I . 



dv 

We obtain the following result. 
Theorem 5.1. We have 



(5.4) 
(5.5) 
(5.6) 



Je[ u e] ~ J[ u o\ 



?1 

2E 



* £ -{Kl + Kh) + 0{e±) 



(5.7) 



as e — > 0. 
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The topological derivative DtJ £ (z) of the potential energy is defined by 

DTMz > := fs,(w)i J ')- (5 - 8) 

where p(e) = ire 2 . So, one can immediately see from (I5.7P that 

D T J £ (z) = -^(K 2 + Kh)- (5.9) 

This formula is in accordance with the one obtained by Novotny et al in [22] (see also |26j). 
In fact, in those papers the plane strain and the plain stress problems are considered, and 
(|5.9p is the formula for the latter problem. 
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A Derivation of ( ED 



The Kelvin matrix (|2.12p can be rewritten as 

d log |x — y 



where 



$ij( x -y) = A'<%log|x-y| +/d'(xi -yi) , i,j = 1,2, (A.l) 

oyj 

A + 3/i , _ A + n 

A — — z — r, 



4vr/x(A + 2/i) n 4tt/z(A + 2 / u) 
Using the operator T(d) defined by (|3.2h one can see that 

d<S> 



dv y 



(x-y) = (r(a y )$(x-y)) J , (A.2) 



or 

(d<5> 

1=1 

We use the formulas 



— (x-y) :=^T, 7 (5 y )^ fc (x-y), k,j = l,2. (A.3) 



<9 , i , Xi-yi 
— log|x - y 1 



"' |x-y| 2 ' 

d 2 , , „ (xi — yi) 2 1 

log x-y = -2 \ yi > + 



I x -y| |x-yr 

<9 2 , (xi — Vi)(xj — Vj) 

— — log x-y = -2 V yi)K 3 U3> iii^j. 
oyiOVj |x - y| 4 
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By (|X3|) . we have 



/ 9$ 

\3Uy J 11 



9 



9 



(A + 2/x)ni- h/un 2 — A'log|x-y| - yi) 

<9yi <9y 2 ' 



d log |x — y| 
dyi 



+ ( ^ n 2^~ + Ani^- ) //(xi - 2/1 ) 
<9yi <9y 2 



9 log |x — y| 



dV2 

\'(\ , o \ <91og|x-y| , w <91og|x-y| 

A (A + 2/i)ni h A /xn 2 i 

oyi oy 2 

, , f , , >, / gjggjx-yj 5 2 log|x-y| 
+ M (A + 2/i)m ^ ^ + (xi - yi) 

, <91og|x-y| , , f s 
n 2 - h n 2 (xi - yi) 



+A/i'(a;i - yi) 



dy2 

d 2 log |x — y| 



dyv 

d 2 log |x — y| 
dyidy 2 



dy 2 



Since A log |x — y| =0 for x / y, we have 



\dUyJ 11 



(A + 2„)(A' - M0n / 1Og i X - y| + MA' - M)- 2 51 ° g ' X " y ' 



dyi 



<9y 2 



, , d 2 log |x — y| 9 2 log |x — y| 

+2/x/x ( ni(xi - yi) hn 2 (xi-yi) 



<9yi 



dyidy 2 



Since 



we obtain 



(A + 2//)(//- AO + 2/i/x' 



Vai/y / 11 



2vr(A + 2//) 
,Oi -yi) 2 



- X') - Afifi 

l x ~y| 



E 

z=i 



- yz 



|x-y 



2 ' 



Similarly, we can compute 



3Vy ) 12 



d 







An 2 — + V n iQ^ ) ( A ' lo g |x- y| + n'(x! - yi) 



Slog |x — y| 



dyi 



( & 9 \ . <91og|x — y| 

+ U 6 "!^ h (A + 2/i)n 2 — /x (xi - yi) 

V oyi dy 2 J dy 2 

X(X' - /i / ) n 2 ^ 1 ° g i X — — + A* (A' - /x')ni^ Y 



9yi 



<9y 2 



/ 9 2 log|x — y| , <9 2 log|x — y| 
+2/i/z ni(xi - yi) ^ — h 2/i/i n 2 (xi - yi) ^-^ 



dy l dy 2 



dy 2 
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Since \(jjf — A') + 2/z// = /i(A' — //), we have 



V <9l/y ) 12 



rw / \A , /] »1 -2/1 , / / ,A ^2 - 2/2 

[A(/x - A ) + Jra 2 - 7* + - A )ni- 



i x - y| 



l x ~y| 



„ / {x 1 -y{f{x 2 -V2) . / (si-yi ){X2-V2) 

"W ra i i m n 2 : 7j 

|x-y| 4 |x-y| 4 

, , {xi-yi){x 2 -y2) ^ 
- 4 w n — 2^ nv 



l x ~y| 



Z=l 



^; - yi 
x - y| 2 



-m(m'-a') 



^2(^1 - 2/1 ) - ni(x2 - 2/2) 



l x -y| 



We also have 



(£) 

V ai/y / 22 



d d \ , <91og|x-y| 

Xn *dy\ +lini W^ {X2 - y2) Wl 











+ Flo h (A + 2/i)ra 2 — A' log |x - y| + //(x 2 - y 2 



dyi 

/<// - A') - 4/x M 



dy 2 

, { X2 - V2) 

|x - y| 2 



d log |x - y| 



E 

z=i 



n zi 



^z - yi 



|x-y 



2 • 



This proves (|3.3p . 
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